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Introduction: We previously modeled a subset of 
domes on Europa with morphologies consistent with 
emplacement by viscous extrusions of cryolava [1]. 
These models assumed instantaneous emplacement of 
a fixed volume of fluid onto the surface, followed by 
relaxation to form domes [1]. However, this approach 
only allowed for the investigation of late-stage eruptive 
processes far from the vent and provided little insight 
into how cryolavas arrived at the surface. Considera- 
tion of dome emplacement as cryolavas erupt at the 
surface is therefore pertinent. A volume flux approach, 
in which lava erupts from the vent at a constant rate, 
was successfully applied to the formation of steep- 
sided volcanic domes on Venus [2]. These domes are 
believed to have formed in the same manner as candi- 
date cryolava domes on Europa [1,3]. In order to gain a 
more complete understanding of the potential for the 
emplacement of Europa domes via extrusive volcan- 
ism, we have applied this new volume flux approach to 
the formation of putative cryovolcanic domes on Euro- 
pa. Assuming as in [1] that europan cryolavas are 
briny, aqueous solutions which may or may not contain 
some ice crystal fraction, we present the results of this 
modeling and explore theories for the formation of 
low-albedo moats that surround some domes (Fig. 1). 
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Figure 1. Candidate cryolava domes with low-albedo moats. 
Domes are 5-6 km wide. 


Approach: An alternative method for inferring 
bulk cryolava rheology, which incorporates the lava 
emplacement stage, is presented here. This method was 
employed in [2] for dome emplacement on Venus. An 
innovative perturbation solution to the generalized 
form of the Boussinesq equation for fluid flow in a 
cylindrical geometry is presented in [4]. The continuity 
equation describing radial expansion of a Newtonian 
fluid with an unbounded (free) upper surface is: 
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[4] found a similarity solution to (1) by transforming 
this equation to an ordinary differential equation when 


a constant volumetric flowrate, Q, at the origin is giv- 
en. This similarity solution is found by perturbation by 
defining a parameter € = 1/(n+1) [2,4]. The new 
independent variable, x, and dependent variable P, that 
transform (1) to an ordinary differential equation are 
defined as: 
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P(x; €) = 4reh*/Q 


where € = % and ® = 1.2 for Europan cryolavas when 
a Newtonian rheology is assumed [2-4]. Following the 
approach from [4], flow thickness, h, as a function of 
time is described by: 


h=h, (PX) (3) 


Results: Fig. 2 shows the solution of a radially 
spreading, Newtonian fluid with a bulk kinematic vis- 
cosity, v = 10° m?/s (equivalent to a bulk dynamic vis- 
cosity, uw = 10° Pas) [1]. Here, the overall “shape” of 
the flow surface and the aspect ratio at t = 40 days is 
very similar to the dimensions of the domes in Fig. 1 
[1, 3]. Note that the bulk kinematic viscosity reported 
above includes the contribution from the brittle cryola- 
va crust [1, 5-6]. The initial viscosity of the erupted 
cryolava may be up to 4 orders of magnitude lower 
than this value [1-2, 5-7], perhaps as low as 10! m7/s. 
In a next iteration of this work, we will explore the 
effects on dome formation while varying v and Q at 
the vent. We will also place bounds on plausible cry- 
olava compositions that would be commensurate with 
the reported viscosity values. 
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Figure 2. Axially symmetric Newtonian fluid flow profiles at five 
times when Q = 550 m/s at the vent 


Moat Formation: Several putative cryolava domes 
are surrounded by low-albedo moats (Fig. 1). A num- 
ber of theories have been set forth for their formation 
Le.g., 3, 8]. Notwithstanding, as these moats are associ- 
ated with putative cryolava domes, it is also possible 


that they were emplaced due to processes associated 
with extrusive cryovolcanism. Cryolava coextrusion 
and lithospheric flexure have thus been explored as 
plausible mechanisms for moat formation. 

Cryolava Coextrusion. In terrestrial volcanic sys- 
tems, magma ascending in dikes and conduits often 
takes the form of two-component flow, in which melt 
of differing viscosities travel side by side to the surface 
[9-12]. As these magmas ascend, viscous segregation 
causes the low-viscosity component to migrate to the 
high-shear regions adjacent to dike walls. This config- 
uration is dynamically stable [9-10, 13-14] and allows 
the low-viscosity component to lubricate the flow of its 
high viscosity counterpart. On Europa, coextrusion is 
likely to take the form of a two component cryomagma 
in which a briny, low-viscosity melt migrates to the 
walls of a fracture, encapsulating a viscous, crystal- 
liquid slurry in the center of the flow where shear is 
negligible. A detailed, analytical treatment of two- 
component effusive flow was undertaken by [12]. Fol- 
lowing this method, the volume flux of cryomagma 
transported during two-component Newtonian flow in 
a cylindrical conduit was derived: 
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G = [g(p; — pr) — AP/H] is the pressure gradient due 
to both the hydrostatic overburden, and the excess 
pressure that would allow melt to be driven from a 
crustal fluid reservoir at depth H [3]. a is the thickness 
of the low-viscosity lubricating layer (i.e., the sheer 
layer) in a conduit of width d. Fig. 3 shows Qwai and 
Qcore aS a function of a for cryomagmas when H = 5 
km and d= 6 m. Here, sdwaii= 20 Pass, Lécore = 4 x 10* Pa 
s and AP = 4 x 10° Pa. For a briny solution ascending 
in an icy crust, p; = 917 kg/m? and p, = 1400 kg/m”. 
For a < d/3, Qcore > Qwali, and the high viscosity melt 
component is rapidly erupted onto the surface (Fig 3). 
This initial eruption stage could therefore provide the 
cryolavas necessary for dome formation at the surface. 
As a approaches d/3, the two-component flow configu- 
ration breaks down and Qvore approaches 0 (Fig. 3). 
The high viscosity cryomagma becomes depleted in 
the conduit, leaving only the low-viscosity component 
available for later eruptions. This low viscosity com- 
ponent may subsequently erupt to form a moat. The 
fluidity of the cryolava involved in this later eruption 
stage would enable it to travel farther from the vent 
than its high-viscosity counterpart, while the moat’s 
low albedo may be a tell-tale sign of its briny composi- 
tion. 
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Figure 3. Volume flux of co-extruding cryomagmas as a function 


of sheer-layer thickness. Here G = 172 Pa/m. 


Lithospheric Flexure. Flexure has been observed 
along fractures on Europa, and ridges, which sit atop 
the elastic lithosphere [15], weigh it down, producing 
secondary uplifts, or forebulges. Galileo imagery re- 
vealed clear cases of crustal warping near ridges, 
which follow the expected flexure geometry for a 
loaded lithosphere. By comparing the observed dis- 
tance of the forebulges with theory, lithospheric thick- 
nesses of 0.2 to 2 km were estimated [16-18]. Extru- 
sion of dome material onto Europa’s surface would 
also place a load on the elastic lithosphere, warping the 
ice shell and causing secondary uplift around domes. 
In the next iteration of this study, dome margins will 
be studied to look for signs of uplift. If a relationship is 
found to exist between candidate cryolava domes and 
lithospheric flexure, it would reinforce a possible cry- 
ovolcanic origin for these domes and suggest that 
moats surrounding these domes are not solely a result 
of material extrusions, but that they may also form 
when the weight of erupted cryolava causes the crust to 
sag into mushy fluid layers beneath. 
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